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Abstract. When the logarithmic slope of the galaxy 
counts is lower than 0.4 (this is the case in all filters at 
large magnitude), the magnification bias due to the lens 
makes the number density of objects decrease, and conse- 
quently, the radial distribution shows a typical depletion 
curve. In this paper, we present simulations of depletion 
curves obtained in different filters for a variety of differ- 
ent lens models - e.g., a singular isothermal sphere, an 
isothermal sphere with a core, a power-law density pro- 
file, a singular isothermal ellipsoid and a NFW profile. 
The different model parameters give rise to different ef- 
fects, and we show how the filters, model parameters and 
redshift distributions of the background populations affect 
the depletion curves. We then compare our simulations to 
deep VLT observations of the cluster MS1008-1224 and 
propose to constrain the mass profile of the cluster as well 
as the ellipticity and the orientation of the mass distribu- 
tion. This application is based solely on deep photometry 
of the field and does not require the measurement of shape 
parameters of faint background galaxies. We finally give 
some possible applications of this method, useful for clus- 
ter lenses. 



Key words: Cosmology - gravitational lensing - galaxy 
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1. Introduction 

Since a few years, a new application of gravitational lens- 
ing in clusters of galaxies has started to be explored, 
namely the depletion effect of number counts of back- 
ground galaxies in cluster centers. This effect results from 
the competition between the gravitational magnification 
that increases the detection of individual objects (at least 
for flux limited data and marginally resolved objects) and 
the deviation of light beam that spatially magnifies the 
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observed area and thus decreases the apparent number 
density of sources. 

This effect was pointed out as a possible a pplica tion of 
the magnification bias by Broadhurst et al. ( |1995 ) where 
they suggested a new method for measuring the projected 
mass distribution of galaxy clusters, based solely on grav- 
itational magnification of background populations by the 
cluster gravitational potential. In addition, they suggested 
that the mass-sheet degeneracy, initially pointed out by 
Schneider and Seitz ( 1995| ) and observed in mass recon- 
struction with weak lensing measures, could be broken by 
using gravitational magnification information which is di- 
rectly provided by the depletion curves. This method has 



been used by Fort et al. J1997D and by Taylor et al. fll998| ) 
to reconstruct a two-dimensional mass map of Abell 1689 
in the innermost 27 arcmin 2 , taking into account the non- 
linear clustering of the background population and shot 
noise. The results are consistent with those inferred from 
weak shear measurements and from strong lensing. How- 
ever, the surface mass density cannot be obtained from 
magnification alone since magnification also depends on 



the shear caused by matter outside the data field (Young 



1981 



Bartelmann and Schneider, 2000). But in practice, 



if the data field is sufficiently large and no mass concen- 
tration lies close to but outside the data field, the mass 
reconstru ction obtained from m agnification can be quite 
accurate ( Schneider et al., 200C ). 

This method is an attractive alternative to weak lens- 
ing because it is only based on galaxy counts and does 
not require the measure of shape parameters of very faint 
galaxies. In addition, it is still valid in the intermediate 
lensing regime, close to the cluster center, without any 
strong modifications of the formalism. Meanwhile, it is 
more sensitive to Poisson noise which increases when the 
number density decreases in the depletion area. Another 



weak point identified by Schneider et al. (2000) is that it 



may significantly depend on the galaxy clustering of back- 
ground objects which can have large fluctuations from one 
cluster to another. 

A second application of the magnification bias that 
has been suggested is to use the shape and the width of 
depletion curves to reveal the redshift distribution of the 
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background populations. This technique was first used by 
Fort et al. ( |1997| ) in the cluster C10024+1654 to study the 
redshift distribution of background sources in the range 
26 < B < 28 and 24 < I < 26.5. They found that 60% ± 
10% of the population in the B band is located between 
z = 0.9 and z — 1.1 while the remaining galaxies are 
broadly distributed around z — 3. The population in the 
I band shows a similar redshift distribution but it extends 
to larger redshift (~ 20% of objects at z > 4). 

The last application of the magnification bias that 
has been explored is the search for constraints on cos- 
mological parameters. This method is based on the fact 
that the ratio of the two extreme radii which delimit the 
depletion area depends on the ratio of the angular dis- 
tances lens-source and observer-source. Consequently it 
depends on the cosmological parameters, as soon as all 
the redshifts are fixed, and mainly on the cosmol ogica l 
constant. This method was first used by Fort et al. ( 1997 ) 
in C10024+1654 and in Abell 370. Their results favor a 
flat cosmology with 0.6 < f^A < 0.85 and are consistent 
with th ose ob tained from other independent methods (see 
White ( 1998 )). This technique requires a good modeling 
of the lens and could be improved by being extended to 
a large number of cluster lenses and by using an inde- 
pendent estimate of the redshift distribution (for example 
from photometric redshifts). However, this method was 
disputed by Asada (1997) who claimed that it is difficult 
to determine £1 by this method without the assumption of 
a spatially flat universe. He also mentioned the uncertainty 
of the lens model as one of the most serious problems in 
this test and concluded that this method cannot be taken 
as a clear cosmological test to determine A. 

In order to better understand the formation of deple- 
tion curves and their dependence with the main character- 
istics of the lenses and the sources, we developed a detailed 
modeling of the curves in various conditions. In Sect. 2 we 
present the lens and counts models used in our simula- 
tions. In Sect. 3 we study the influence of the lens mass 
profile, with several sets of analytical mass distributions. 
In Sect. 4 we explore the influence of background sources 
distribution, selected through different filters. An appli- 
cation on real data is presented in Sect. 5, in the cluster 
MS 1008-1224, observed with the VLT. We compare these 
observations with the simulations in order to derive some 



constr; hits on the mass profile and the cllipticity of the 
cluster, as well as on the background sources distribution. 
Some prospects about the future of this method are given 
as a conclusion in Sect. 6. 

Throughout the paper, we adopt a Hubble constant of 
H = 50 km s" 1 Mpc -1 , with A = and Q = 1- 



2. Modeling depletion curves 



2.1. The magnification bias 

The projected number density of objects magnified by a 
factor /i through a lensing cluster and with magnitude 
smaller than to can be written in the standard form: 

N(< to) = N (< to) fi 2 -^'- 1 

where No is the number density of objects in an empty 
field and 7 is the logarithmic slope of the galaxy num- 
ber counts. Note first that this relation applies only for 
background galaxies. This means that we implicitly work 
at faint magnitude where the foreground counts are re- 
duced and do not affect the slope 7. Note also that the 
magnification // depends on the source redshift through 
the geometric factor T)ls/^os m a significant way. For 
example, for zl = 0.4, Dls/^os changes by more than 
40% when z$ varies from 0.8 to 2. For a lower redshift 
cluster lens, the effect should be reduced, at least for all 
sources at z > 0.8. 

An exact writing of the magnification bias, which also 
takes into account the local effect of magnification, chang- 
ing with radius r within the lens can be re-formulated as: 

N(r, < to) = 

N (< to) 1 ' 

f* L n(z, <m)dz + /+°° n(z, < to) ^(r, z) 2 - 5 ^ 1 dz 
Jo + °° n(z, < to) dz 
r+°° n(z, < to) [1 - n(r, z) 2 -^- 1 } dz 
N (< to) 

where n(z, < to) is the density of galaxies of apparent 
magnitude smaller than to, located at redshift z. "f(z) rep- 
resents the shape of the luminosity function at redshift z 
and we will assume in the following that 7 does not de- 
pends on redshift, i.e. most of the faint galaxies are repre- 
sentative of the power-law part of the luminosity function. 
In addition we suppose that the slope is a constant with 
redshift. 

It is clear that the radial behavior of the ratio N(r, < 
to)/No(< to) is strongly related to n(r). It decreases from 
the center, up to a minimum at r ~ Re when magnifica- 
tion goes to infinity (Einstein radius) and then increases 
up to 1 again, when magnification goes to 0. This is the 
so-called radial depletion curve which has been detected 



1998; Athreya et al., 2000) 



in a few clusters already (Fort et al., 1997; Taylor et al 



From the lensing point of view, the magnification \i 
can be written as (Schneider et al., 1992) : 

1 1 

A* = -7— T = T, ..2 (2) 



det.4 (1 - k) 2 



'7 



where A is the magnification matrix. The convergence k 
and the shear 7 are expressed in cartesian coordinates as 
a function of the reduced gravitational potential ip by : 

_ 1 /9V #V 
K ~ 2 \8x 2 + dy 2 



(3) 
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7 = V 7i + ll 



(4) 



with 



7i 



1 ( Ftp _ av 

2 \<9x 2 9y 2 



and 



72 



d 2 <p 

dxdy 



All these factors depend on the lens and source redshifts, 
as <f is related to the true projected potential $ by : 



2 1 Dr.,D. 



C 2 Dos 

where £ is a characteristic scale of the lens. 

In the weak lensing regime we have the approximation 
ji{r) ~ l+2/t(r) flBroadhurst et al., 1995). So the potential 



interest of the depiction curves is that they trace directly 
the convergence distribution, or equivalently the surface 
mass density distribution. In principle, they may allow 
an easy mass reconstruction of lenses in the weak regime. 
We will see below what are the main limitations of such 



curves, already explored by Athreya et al. (200C) 



2.2. The cluster lens 

Different sets of models are used, with increasing complex- 
ity. In all cases we suppose that the lens redshift is 0.4, 
with reference to the cluster C10024+1654 studied by Fort 
et al. ( |1997| ) . With the adopted cosmology this means that 
the scaling corresponds to 6.4 kpc for 1". For each 
model we will use an analytic expression and develop it 
to compute analytically the magnification and its depen- 
dence with redshift and radius. All the models are scaled 
in terms of the Einstein radius Re- 

2.2.1. The singular isothermal sphere (SIS) 

This model is the simplest one which can describe a clus- 
ter of galaxies and is very useful for the analytical calcu- 
lations of the gravitational magnification. However it has 
also some physical meanings: 

• it corresponds to a solution of the Jeans equation and 
thus it can be written as a function of the observed 
velocity dispersion; 

• in a non-collisionnal description of the collapse of a 
self-gravitating system, the violent relaxation, dur- 
ing which particles exchange energy with the average 



field, leads to an isothermal distribution (Binncy and 



Tremaine, 1987) 



But, this model is only valid inside certain limits due to 
the divergence of the central mass density and of the total 
mass to infinity. This divergence has no consequences on 
our work because the validity limits of the model corre- 
spond to the limits of the depletion regime. The density 
of matter can be written as : 



p(r) = Pe 



r 

Re~ 



(5) 



where Re is the Einstein radius : 

o 2 D ls Dql 
Re = 4?r — — 

c 2 D s 

and a is the velocity dispersion along the line of sight. The 
gravitational magnification is simply : 

Re 



^ = 1 - 



(6) 



2.2.2. The isothermal sphere with core radius 



This model avoids the divergence of the mass density in 
the inner part of the cluster with an internal cut-off of the 
density distribution. We chose the following distribution, 
as described in Hinshaw and Krauss (1987) or Grossman 
and Saha (1994). The density of matter is given by : 



P( T ) = 7T- 



1 



2tt GR? c (1 



A 2 



1 



|A 2 -1|)^ 



(7) 



RcVZ 1 



1, 



where Rc is the core radius, x — r/Rs and $o is the 
central value of the gravitational potential 
In this case, the Einstein radius is Re 
where : 

_ Aira 2 D LS Dql 
c 2 Dos Rc 
The magnification can be written as : 



1 - A (l + x 2 | A 2 - l| 



(8) 



1- A 1 



1 



Physically, in most cases the core radius is smaller than 
or comparable to the Einstein radius, so we do not expect 
strong effects in the outer parts of the depletion curves. 

2.2.3. A power-law density profile 

This model is a generalization of the SIS. It can be used 
in order to test the departure from an isothermal profile 
far away from the cluster center. The density is given by : 



p = PE 



Re 



(9) 



where a is the logarithmic slope of the density profile. In 
order to keep a physical model for the mass distribution, 
we must have a < 3. The Einstein radius can be written 



Re = \ -o 



4 D LS D OL 



Dos 

4np E R E / (3 



GMe I a /2 



where M E — 4:Trp E R E / (3 — a) is the integrated mass in- 
side the Einstein radius and Jg = J Q °°(1 + ^ 2 )~' 9 c?£. The 
magnification is then : 



P~ X = 



l-(2-a) 



Re J 



1 {re) 



(10) 
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2.2.4. The singular isothermal ellipsoid 



This t^i pe of model is mterestigg as a lot of clusters have 



elliptic al shapes in their galaxy distribution or their X- 



ray isophotes ( 


Buotc and Canizares, 1992 




Buotc and 


Canizares, 1996 


; Lewis et al, 1999; |Soucail ct al, 2000]). 



Thus, an elliptical lens represents a more realistic model 
although it is still reasonably simple. In order to study the 
effect of the ellipticity of the potential e = 1 — - (a > b) 
on the depletion, we introduced the singular isothermal 
ellipsoid in our simulations ( Kormann et al., 1994 ). Using 
polar coordinates (r, ip) in the lens plane we introduce 



( = \J (rcostp) 2 + (1 — e) 2 (rsinipY 



(11) 



which is constant on ellipses with minor axis (, and ma- 
jor axis C/ (1 — e)- The surface mass density can then be 
written as : 



£ (r, if) 



VT^la 2 1 
2G C 



(12) 



where a is the velocity dispersion along the line of sight. 
The convergence is : 



. = ^!^(^)x 4 - G DolDls 



(13) 



^critical ' C 2 DoS 

The magnification writes then quite simply: 

fi~ 1 = l-2K (14) 

2.2.5. NFW profile 

We introduced in our simul ations the universal density 
profile of Navarro et al. ( 1996 ) for dark matter halos which 
was found from cosmological simulations of the growth of 
massive structures. This profile is a very good descrip- 
tion of the radial mass distribution inside the virial radius 



T2oo- Wright and Brainerd (200C) have compared it with 
a SIS for several cosmological models. They find that the 
assumption of an isothermal sphere potential results in 
an overestimate of the halos mass which increases linearly 
with the value of the NFW concentration parameter. This 
overestimate depends upon the cosmology and is smaller 
for rich clusters than for galaxy-sized halos. 
The NFW density profile is given by : 



p(r) 



0~cPc 



(r/r s ) (1 + r/r s y 



where r s = r^jc is a characteristic radius and p c is the 
critical density, c is the concentration parameter and 5 C is 
a characteristic over-density for the halo which is related 
to c by the requirement that the mean density within the 
virial radius should be 200p c . This leads to : 



200 



3 ln(l + c) - c/(l + c) 



If we take x — r/r s , the surface mass density and the shear 



of a NFW lens can be written as (Wright and Brainerd 



2000) 



knfw(z) 



7nfw(^) 



r s $cPc 



r s S c p c 



f<(x) (* < 1) 

2/3 (x = 1) 

f>(x) (x>l) 

g<{%) (x < i) 

10/3 + 41n(l/2) (x = l) 

9>(x) (x > 1) 



where E c is the critical surface mass density. /< (x) and 
f>(x) express as : 



f<(x) = 
f>(x) = 




(x 2 - 1) 

while g < (x) and g > (x) express as 
!ar S th \/S 4 

9<{x) 




(x 2 -l) (x 2 -l)yT^ 



9>(x) 



8 arctan 




The magnification by the NFW lens is then : 

-1 /1 \2 2 

P = (1 — «nfwJ — 7nfw 



(16) 



2.3. The galaxy redshift distribution 
2.3.1. Analytical distribution 

We used the analytical redshift distribution introduced by 



Taylor et al. (1998) to fit the redshift distribution of the 
galaxies in the range 20 < / < 24 for 0.25 < z < 1.5 : 



(15) N(z,I< 24) = N {I< 24) 



*2r(3/a; 



exp 



(17) 



with a = 1.8 and z* = 0.78. Here, we extend the redshift 
range between z = 0.05 and z = 3 and we take Nq(I < 
24) = 12 objects per arcmin 2 which is the count rate found 



by Taylor et al. ( |1998[ ) for red galaxies with 20 < / < 24. 
The distribution shows a maximum at z = 0.8 and has 
an average redshift of < z >~ 1. It is worth noting that 
with these parameters, the redshift distribution does not 
correspond to very deep observations. 



Christophe Mayen & Genevieve Soucail: Depletion curves of galaxy number counts 



5 



Table 1. Magnitude ranges adopted for the field num- 
ber counts in the simulations (see text for more details 
concerning the number count model). 



Filter 


Magnitude range 


Faint end slope 


B 


22 - 28 


0.14 


I 


20 - 26 


0.23 


K 


18 - 23 


0.22 



2.3.2. The magnitude-redshift distributions in the U to K 
photometric bands 



We also used a model of galaxy number counts from 
which the magnitude-redshift distribution in empty field 
can be computed with a large set of photometric bands. 
This model was developed by Bezecourt et al. ( jl998 ) and 
largely inspired by Pozetti et al. (1996). It includes the 
mod el of g alaxy evolution developed by Bruzual and Char- 
lot ( 1993 ), with the upgraded version so-called GISSEL, 
and standard parameters for the initial mass function 
(IMF) and star formation rates (SFR) for different galaxy 
types. In order to reproduce the deep number counts of 
galaxies, evolution of the number density of galaxies is in- 
cluded to compensate for the smaller volumes in an f2o = 1 
universe, followi ng the prescriptions of Rocca-Volmerange 
and Guiderdoni ( 199C ) . This model reproduces fairly well 
deep number counts up to B = 27 or / = 25, as well as 
r cdshift distributions o f galaxies up to B = 24 or / = 22 
( Bezecourt et al., 1998 ). 

The advantage of this model is that the galaxy distri- 
bution can be computed for any photometric band, and 
the effects of the color distribution of sources can be ex- 
plored. In practice we will limit our study in this paper to 
B, I and K bands (Table Q). 



3. Influence of the lens parameters on depletion 



For our first set of simulations, we used the analytical red- 
shift distribution for the sources, allowing a fully analyt- 
ical treatment of the simulations. The contamination by 
foreground objects with this distribution is weak (about 
~ 9% of foreground galaxies for zl = 0.4). In each case, 
we computed depletion curves, and then examined their 
behavior. Note that we limited our analysis outside the 
first 20" from the center where the signal cannot be con- 
strained observationally (decrease in the observed area 
for each point, obscuration by the brightest cluster galax- 
ies . . . ). For these reasons, this "forbidden" area will be 
shaded in each plot. 




a = 


1800 


km 


s"' 


o = 


1600 


km 


s -1 
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1400 


km 


s _1 




1 200 


km 
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(7 = 


1000 


km 


s" 1 
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Radial distance from the cluster center (arcsec) 



1400 
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Fig. 1. Top: Depletion curves obtained for different ve- 
locity dispersions ranging from 1000 to 1800 km s . Bot- 
tom: Intensity of the minimum (solid line), position of 
the minimum in Mpc (dotted line) and half width at half 
minimum in Mpc (dashed line) of the depletion curves as 
a function of the velocity dispersion. 



3.1. Influence of the velocity dispersion 

We used a SIS model and varied the velocity dispersion 
a from 1000 to 1800 km s" 1 (Fig. 0). As expected the 
radius of the minimum increases with a. Indeed for a SIS, 
the Einstein radius scales exactly as er 2 , so the minimum 
of the depletion curve, which roughly corresponds to the 
maximum Einstein radius (at z ~ 4 or more), also scales 
as a 2 . 

More surprisingly, the depth of the depletion at the 
position of the minimum is roughly constant and, at first 
order, does not depend on the velocity dispersion of the 
cluster or its total mass. Part of this effect is due to the 
density of foreground galaxies, but part of it is intrinsic 
to this SIS model, as other potential shapes do not show 
this property (see below) . On the contrary there is a clear 
dependence on the half width at half minimum with a 
(Figure0). 
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Fig. 2. Depletion curves obtained with a = 1400 km s _1 
for different core radii ranging from 25 to 125 h$Q kpc. 
The value of the Einstein radius at the cluster redshift is 



200 h7} kpc 



3. 2. Influence of the core radius 

The introduction of a core radius in the model (Fig. ||) 
does not affect significantly the outer region of the de- 
pletion area and has little effect on the position of the 
minimum. As Rc increases, the inner width of the curve 
is enlarged and its slope is decreased. In fact, the study 
of this area for our purpose has no interest, because the 
spatial resolution of galaxy number counts variations is by 
far larger than this scale. In addition, in rich clusters of 
galaxies, the central density of galaxies is large enough so 
that empty areas are quite small between the envelopes of 
large and bright galaxies. 

In practice, the only observable effect of a core radius 
is a deformation of the inner depletion curve with a sig- 
nificant departure from a symmetry with respect to the 
outer part. But we suspect that a quantitative estimate of 
Rc would be difficult to extract from the signal. 

3. 3. Influence of the slope of the mass profile 

The simulations are done with the power-law density pro- 
file, characterised by its slope a (Fig. |^). The minimum 
position of the depletion area does not depend significantly 
of a, because our scaling of the potential was fixed at 
the Einstein radius, nearly independent of a. On the con- 
trary, the two other typical features of the depletion area 
(half width at half minimum and intensity of the min- 
imum) strongly depend on this parameter. The width of 
the curve represents roughly the mass dependence with ra- 
dius. In the case a = 1.7 for example, the shallower slope 
of the mass radial dependence creates a larger width of 
the depletion curve. The half width at half minimum of 
the curve decreases when a increases. The reason is the 
same as with the SIS, that is to say, when a increases, the 




Fig. 3. Top: Depletion curves obtained with pe = 2.25 x 
10 15 M0/Mpc 3 and for different slopes a ranging from 1.7 
to 2.1. Bottom: Intensity of the minimum (solid line), 
position of the minimum in Mpc (dotted line) and half 
width at half minimum in Mpc (dashed line) of the deple- 
tion curves as a function of the slope a. 



mass is concentrated in the inner part of the cluster and 
the depletion effect is less extended to the outer regions. 

3.4. Influence of the ellipticity of the potential 

In these simulations, we use the singular isothermal el- 
lipsoid potential, with a = 1400 km s _1 , and we study 
the depletion curves along the minor and major axis (Fig. 
^) . The increase of e does not affect the minimum value of 
the depletion area but leads to an increase/decrease of the 
half width at half minimum and of the minimum position 
along the major/minor axis as an homothetic transfor- 
mation. The relative positions of the two minima of the 
depletion area along the main axis gives immediately the 
axis ratio and thus the ellipticity of the potential. An ap- 
plication of this differential effect along the two main axis 
is presented in Sect. |5[ 
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Fig. 4. Top: Variation of the depletion curve with el- Fig. 5. Top: Depiction curves obtained with c = 15 and 
lipticity along the minor axis. Bottom: Idem along the for different virial radii. Bottom: Intensity of the mini- 
major axis. The streching of the curves with ellipticity is mum (solid line), position of the minimum in Mpc (dotted 



characterised by an nomothetic transformation. 



line) and half width at half minimum in Mpc (dashed line) 
of the depletion curves as a function of the virial radius. 



3.5. Study of NFW profile 

With the NFW model we varied separately the virial ra- 
dius and the concentration parameter. We have chosen the 
values for these two parameters in the range of those found 
by Navarro et al. (1996). For the first set of simulations, 
we adopted c = 15 and varied r2oo from 1800 to 3600 
kpc (Fig. ^). An increase of the virial radius (more mas- 
sive cluster) leads to an increase of the three characteristic 
features of the depletion area and to the appearance of a 
bump in the central region which grows with r2oo- The 
increase of the half width at half minimum of the deple- 
tion curves with r2oo can be explained by the fact that 
when r2oo increases there is more mass in the outer re- 
gions. Consequently, the depletion area is more extended 
to the outer part of the cluster and reaches its asymptotic 
value less rapidly. The second set of simulations was done 
with r2oo = 3000 kpc and we varied c from 7 to 25 (Fig. 
^). The variation of c affects only the inner part of the 
depletion curves. An increase of the concentration param- 



eter (cluster with smaller characteristic radius) leads to 
an increase of the intensity and position of the minimum 
with a steeper slope for the small values of c. Contrary 
to r2oo, the variation of c does not affect significantly the 
half width at half minimum of the depletion area. 

4. Influence of background galaxies on depletion 
curves 

In this section, we use a unique model for the lens, i.e. a 
SIS profile with a velocity dispersion of 1500 km s _1 for 
a cluster located at zl = 0.4, to avoid effects due to the 
lens. 

4-1- Effect of the lens on the magnitude- redshift 
distribution of background populations 

Fig. shows the distortion of the magnitude-redshift dis- 
tribution in the B band along the curve. The results are 
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Fig. 6. Top: Depletion curves obtained with T2oo = 3000 
kpc and for different concentration parameters. Bottom: 
Intensity of the minimum (solid line), position of the mini- 
mum in Mpc (dotted line) and half width at half minimum 
in Mpc (dashed line) of the depletion curves as a function 
of the concentration parameter. 



of the depletion area (Fig. ||). The magnitude limits of 
the simulations correspond to a typical observation time 
of about 2 hours on an 8-meter telescope with good imag- 
ing facilities. The contamination by foreground objects 
(z < 0.4) in the 3 filters is weak (~ 8% in the I and 
K bands and ~ 6% in the B band for which the counts 
are deeper). 

The evolution of the position of the minimum of the de- 
pletion area reflects the evolution of the mean redshift of 
the different populations (< z >~ 1.57 for the K band, 
< z >~ 1.76 for the B band and < z >~ 1.79 for the 
I band), but with a weak effect. The variation of the 
intensity of the minimum of the depletion area is con- 
nected for one part to the fraction of foreground objects 
and for the other part to the slope 7 of the luminosity 
function. As this slope decreases from I band to B band 
(7/ = 0.23, j K = 0.22,7s = 0.14) with the magnitude 
thresholds we have chosen, the intensity of the minimum 
of the depletion area decreases in the same way from I 
band (~ 0.31) to B band (~ 0.26). Both features of the de- 
pletion curves are in fact weakly dependent on the choice 
of the filter and may be difficult to distinguish with real 
data, at least while the variation of the redshift distribu- 
tion with wavelength is rather smooth. 
On the contrary, the half width at half minimum of the 
depletion curve is affected by the redshift distribution of 
background populations in a more sensitive way. The in- 
crease of the half width at half minimum from B band 
to K band is anti-correlated with the fraction of objects 
at high redshift (z > 2) for these populations which also 
decreases in the same way from B band (~ 37%) to K 
band (~ 29%) and seems to reflect essentially the bulk of 
galaxies at redshift around 1. 



5. Depletion effect in the cluster MS1008 1224 



qualitatively the same in the others bands. Near the clus- 
ter center, where the effects of the gravitational magnifica- 
tion are stronger, objects located just behind the lens show 
a gain in magnitude larger than 1 . In the minimum area of 
the curve, the competition between gravitational magnifi- 
cation and dilatation is more favorable to high z objects 
(z > 2.4). In the outer part of the cluster, the effects of the 
gravitational magnification decrease and one finds again 
the distribution in empty field (asymptotic limit of the 
depletion curves). 

4-2. Differential effects in several filters 

We first studied the evolution of the depletion area with 
the redshift distribution of the background population by 
considering several filters. We tried to determine some 
characteristics of the redshift distribution of background 
population which could be connected to typical features 



We have applied our simulations to the cluster MS1008- 
1224 (z = 0.3062, Lewis et al. Ql999| )) which is one cluster 
of the E instein Observatory Exten ded Medium Sensitivity 
Survey ( Gioia and Luppino, 1994 ) . It is a very rich galaxy 
cluster, slightly extended in X-rays, and it is also part of 
the Canadian Network for Observational Cosmology Sur- 
vey ( Carlberg et al., 1996 ). Its galaxy distribution is quite 
circular, surrounding a North-South elongated core. There 
is a secondary clump of galaxies to the North. The X-ray 



luminosity is L x (0 .3-3.5 keV)= 4.5 x 10 44 erg.s" 1 flGioia 



|and Luppino, 1994 ), the X-ray temperature inferred from 
ASCA obser vations is Tx = 7.3 keV ( Mushotzky and 
[gcharf, 1997) and the radio flux at 6 cm is lower than 
0.8 mJy (Gioia and Luppino, 1994). Some gravitation- 
ally lensed arcs to the North and to the East of the field 



have been reported by Le Fevre et al. (1994) as well as by 
Athreya et al. d2000[). 
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Fig. 7. Upper-left: Magnitude-redshift distribution of the galaxies in empty field computed from the model of 
Bezecourt et al. (1998), in the B-band; Upper-right: Depletion curve for rnBUm = 28. The open squares show 
the points where the following Ni ens (m, z,r) distributions are computed. Middle-left: r = 31" (1); Middle-right: 
r = 42" (2); Lower-left: r = 70" (3); Lower-right: r = 125" (4). 



5.1. The observations 

Data were obtained during very deep observations 
with FORS and ISAAC during the science verifica- 



t ion phase of the VLT-ANTU (UT1) at Cerro Paranal 
( http : //www, eso . org/ science/utlsv ). Multicolor pho- 
tometry was obtained in the B, V, R and I bands with 
FORS1 (field-of-view 6.8' x 6.8') and in the J and Ks bands 
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Fig. 8. Depletion curves obtained in different filters for a SIS profile with a = 1500 km s _1 (top left) and corresponding 
field redshift distributions for each filter and magnitude limit. The shaded area on the redshift distributions indicates 
the foreground contaminating population. 



with ISAAC (field-of-view 2.5' x 2.5') with sub-arcsecond 
seeing in all cases. The total exposure times are respec- 
tively 2880 seconds in J, 3600 seconds in Ks, 4050 seconds 
in I, 4950 seconds in B, 5400 seconds in V and R. We 



used the SExtractor software (Bertin and Arnouts, 1996) 
to construct our photometric catalogue. The completeness 
magnitudes are respectively J = 24, Ks = 22, I = 25.5, 
V = 26.5, B = 26.5 and R = 26. Our values are identical 
to those of Athreya et al. ( |2000 ) concerning the FORS ob- 
servations but they are one magnitude deeper concerning 
the ISAAC ones. In practice, we used only the FORS data 
which extend to a larger distance and are more suited to 
our study. 

In order to remove part of the contamination by cluster 
members, we identified the elliptical galaxies by their po- 
sition on the color-magnitude diagram (Fig. |l0j). So back- 
ground sources were selected by removing this sequence, 
essentially valid at relatively bright magnitudes. This sta- 
tistical correction is not fully reliable as it does not elimi- 



nate bluer cluster members or foreground sources. But it 
partly removes one cause of contamination, especially sig- 
nificant close to the cluster center. Our counts may also 
present an over-density of objects in the inner part of the 
cluster (r < 20") due to the non-correction of the surface 
of cluster elliptical galaxies we have removed, effect which 
is more sensitive in the inner part of the cluster where 
these galaxies are dominant. But as it is not in the most 
interesting region of the depletion area, we did not try to 
improve the measures there. 

In addition, due to the presence of two 11 magnitude 
stars in the northern part of the FORS field, two occulting 
masks were put to avoid excessive bleeding and scattered 
light. We took into account this partial occultation of the 
observed surface for the radial counts above a distance of 
130" (see Fig. ^|). This may possibly induce some addi- 
tional errors in the last two points of the curves, which 
are probably underestimated, because of the difficulty to 
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the barycenter of the points with the lower density around 
the cluster. We did not attempt to fit it and preferred to 
fix it 15" North of the cD, following both the X-ray center 



Fig. 9. VLT FORS1 image of MS1008-1224in the B band. 
North is up, East is left. The field-of-view is 6.8' x 6.8'. 
Note the two occulting masks hiding the brightest stars in 
the field. 




Fig. 10. Color-magnitude diagram R-I versus R of 
MS1008-1224. The two horizontal lines delimit the se- 
quence of the cluster elliptical galaxies, which has been 
removed from our catalogue. 



position (Lewis et al., 1999) or the weak lensing center 
( [Athrcya et al., 200C| ). 

The radial counts were performed in a range of 3 mag- 
nitudes up to the completeness magnitude in the B, V, R 
and I bands and up to a radial distance of 210" from the 
center, which covers the entire FORS field. The count step 
was fixed to 30 pixels (6") in the innermost 80" and for 
the rest of the field we adopted a count step of 60 pixels 
(12") to reduce the statistical error bars. These values are 
a good balance between statistical errors in each bin which 
increase for small steps and a reasonable spatial resolution 
in the radial curve, limited by the bin size. The depletion 
curves obtained in the B, V, R and I bands are shown on 
Fig. [nj At some small limiting radius (r < 20"), the ratio 
between the area of the count rings and the area covered 
by the cluster galaxies is close to unity. For this reason, 
we exclude the measures obtained in these innermost rings 
for any study of depletion effects (shaded regions on figure 
0). 

We fitted the observed depletion curve in the I-band 
with three of the mass models discussed above, namely a 
singular isothermal sphere, a power-law density profile and 
a N FW p rofile, and with the galaxy distribution described 
in § 2.3. 2\ For each model, a \ 2 minimization was intro- 
duced to derive the best fits and their related parameters 
(Table ||). The first fit included all the data points (x|p) 
and gave poor results with a reduced \ 2 always close to or 
larger than 2. A second fit was done after removing some 
clear deviant points (x^a)'- the last two points probably 
poorly corrected from edge effects, and those associated 
to the overdensity seen at r ~ 80". This bump is easily 
identifiable in the V, R and I curves, and can be partly 
explained by the presence of a background cluster lensc d 
by MS1008-1224 and identified by Athreya et al. fl2000D . 
Nevertheless, even if we remove from our data all the lower 
right quadrant of the field where this structure is located, 
the bump is still there although significantly reduced. This 
suggests that it may be more extended behind the clus- 
ter center than initially suspected. This second fit gave 
more satisfying results. So we will consider the parame- 
ters associated with the best fit of each model as our best 
results. 



estimate the surface of the masks and some edge effects 
at the limit of the field. 



5.2. Depletion curves and mass density profile 

The first step would be to locate the cluster center. Its 
position is rather difficult to estimate directly from the 
distribution of the number density of background sources, 
although in principle one should be able to identify it as 



The velocity dispersion derived from the SIS model 
km s _1 ) is in good agreement with 



1200 



+200 
-175 



Ht 

the value measured by Carlberg et al. ( 1996| ) (ct ^ s = 



1054 ± 107 km s ) but is more discordant with the 
value of 900 km s _1 inferred from the shear analysis of 
Athreya et al. ( |200C| ). 

The slope of the potential fit with a power-law den- 
sity profile is close to an isothermal one (a = 1.88), 
although slightly shallower. 
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Fig. 11. Depletion curves measured in MS1008-1224 in several spectral bands up to the completeness magnitude of 
the catalogues. For each plot, the best fit by the three models described in text are given. Error bars correspond to 
Poisson statistical noise. The shaded area on each plot corresponds to regions where the study of the depletion does 
not make sense (see text for more details). 



• For a NFW profile, we find a virial radius (Y200 = 
3.2/i^q 1 Mpc) and a concentration parameter (c = 8.9) 
quite in good agreement with those of Athreya et al. 

derived from weak lensing measures. 

• The comparison between the 3 fits favors a NFW pro- 
file as the best fit of our depletion curves (Table |), also 
in agreement with the shear results for this cluster. 

Note however that we did not use the B-band data in 
the fit because the depletion effect is less obvious and much 
noisier than in the 3 other filters. This is due to the fact 
that for this magnitude range the logarithmic slope of the 
counts is close to the critical value 0.4 and the depletion 
effect is strongly attenuated. Probing fainter objects in 
this band (26 < B < 28), even beyond the completeness 
limit, strengthens the depletion signal again , as it is shown 
in Figure [l2| In addition this deeper magnitude range may 
be used to scan a different redshift distribution with a 
higher density of high redshift objects. This is one of the 



future prospects of the depletion curve analysis to perform 
a multi-wavelength analysis to try to disentangle lensing 
effects from the statistical distribution of the sources. 

5.3. Ellipticity and orientation of the mass distribution 

We have also used the depletion effect to determine the 
orientation of the main axis and the ellipticity of MS1008- 
1224. For this purpose, we divided the area of the FORS 
field into four quadrants where quadrants 1 and 2 are in 
opposition, and we performed the counts for radial dis- 
tances from the cluster center to 130" (Fig. The posi- 
tions of the four quadrants are marked by 0, the clockwise 
angle between the North-South axis and the median of 
quadrant 1. We then computed the depletion curves cor- 
responding to the counts in quadrants 1+2 and 3+4. We 
varied 9 and followed the shift of the minimum position of 
the depletion area of the two curves. We can measure the 
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Table 2. Results of the fit of the model parameters with 
a x 2 minimization. The errors are given at the 99.9 % 
confidence level. 



Model 



Parameter 



reduced x 



SIS 
PLP 
PLP 
NFW 
NFW 



SIS 
PLP 
PLP 
NFW 
NFW 



All points 

a = 1200±^° km s _1 
p E = 2.4 ±8;| xlO 15 M /Mpc 3 
a = 1.88±8;?| 

r 200 = 3.2±8l4 Mpc 
c = 8.9±|;g 

Without 3 deviant points 

cr = 1200±i?g km s" 1 
p B = 2.4±Jjx 10 15 M /Mpc 3 



T2U0 



= 1.88±8:as 
3.2±8:| Mpc 

= 8.9±!:S 



1.85 
1.87 
1.89 
2.16 
2.18 



0.49 
0.62 
0.65 
0.44 
0.45 



orientation of the two main axis when the shift between 
the two minima is maximum. For this value of 0, the rel- 
ative positions of the two minima give the main axis ra- 
tio and consequently the ellipticity. The results, obtained 
from I-band data only, are similar in the others bands. 
From the I-band data only, we found 8 = —10° ± 3° and 
an ellipticity of e = 0.18 ± 0.04 (Fig. ||). These results 
agree well with the mass map orientation derived from 
weak lensing as well as the galaxy number density map 
( Athreya et al., 200C| ). A small discrepancy arrises with 
the X-ray gas distribution which orientation is shifted a 



few 10 to 20 degrees West from our measure ( Lewis et al 
1999|)~|But globally, all these results are consistent within 
each other. 




50 100 150 

Radial distance from the cluster center (arcsec) 




Fig. 12. Top: Depletion curve obtained for 26 < B < 28. 

The fits by the three models described in text are given. 
Bottom: Redshift distribution of background objects 
which corresponds to the best fit of the above depletion 
curve. 



6. Conclusions 

The gravitational magnification by cluster lenses repre- 
sents an effective tool to probe the distant universe and the 
mass distribution in the lenses. The first aim of this paper 
was to study the effects of lens mass profile model param- 
eters on the typical features of depletion curves. We have 
also attempted to characterize some features associated 
with the background redshift distribution of the galaxies. 
Models were constructed in three bands covering a large 
spectral range and by using five different lens models. 

Our simulations agree well with very deep and high 
quality images of the cluster MS1008-1224 obtained with 
the VLT and FORS. The depletion effect is clearly seen 
in this cluster, and we have fitted its radial variation with 
several sets of mass profiles. The results are quite sat- 
isfying as we are able to constrain the mass profile up 
to a reasonable distance from the center (about 200", or 
equivalently 1.1 h^Q Mpc). Our results marginally favor 



the NFW mass profile over the isothermal profile, and 
more significantly reject a power-law distribution, essen- 
tially because of the steep rise of the depletion curve just 
after its minimum. Note that this region corresponds to 
the "intermediate" lensing regime, where shear measure- 
ments are more difficult to relate to the mass distribu- 
tion. We have also studied the shape of the depletion area, 
which is easy to relate to the ellipticity of the mass distri- 
bution. We were then able to constrain the ellipticity and 
orientation of the potential of MS1008-1224 with a good 
accuracy. 

This preliminary study highlights the need for addi- 
tional exploration of several issues not fully explored in 
the present paper. For example, the question of clustering 
of the background sources still remains, although in the 
case of MS1008-1224, we have shown how it can be partly 
eliminated. Schneider et al. (2000) also mention this prob- 
lem and insist on the fact that, at deep magnitudes, the 
two point correlation function of the sources is still quite 
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Fig. 13. Depletion curves obtained in the I band for 6 = 
0° (top), 6 = -10° (middle) and 6 = -20° (bottom). The 
radial step has been reduced to 6" only for a more accurate 
localisation of the minima 



tion does not exceed a few percent for an angular separa- 
tion of 10", typical of our bin size. Previous measurements 
were less optimistic in the sense that their values attained 
8 to 10 % at 10", even at very faint magnitudes. But al- 
though there is still only poor knowledge of the amplitude 
of the correlation function at the faintest levels, we can 
hope that the effect will not be dramatic in our studies 
of the magnification bias, provided we retain a reasonable 
bin size at several arcseconds, to wash out most of the 
inhomogeneities . 

Following prelimi nary approaches from the observa- 
tiona l point of view ( [Taylor et al., 1998 ; Athreya et al 



2000 ) or a more theoretical one ( Schneider et al., 2000 ) 



one now clearly needs to extensively and quantitatively 
compare the weak lensing approach with the magnifica- 
tion bias. In particular, with the new facilities of deep 
wide-field imaging presently available, most of the difficul- 
ties related to a small field-of-view can be overcome: for 
weak lensing measurement, a complete mass reconstruc- 
tion requires shear measurements up to the "no-shear" re- 
gion in the outer parts of the cluster to integrate the mass 
inwards. The absolute normalization of the field number 
counts for the magnification bias can also be estimated 
outside the cluster, in exactly the same observing condi- 
tions (filter, magnitude limit, seeing, . . . ), giving an abso- 
lute calibration of the depletion effect. 

Moreover, the full 2D mass reconstruction of the clus- 
ter from the depletion signal alone should be tested. As 
the signal is directly related to the magnification fi(r), and 
then to k(t) and j(r), it should be in principle feasible 
to invert the depletion map to produce a non-parametric 
mass map. In practice, the reconstruction is simpler as 
soon as the shear 7 becomes small, because in that case 
ll and k are simply linearly related. In any case, we have 
shown in this paper that it is rather easy to reach the outer 
parts of the cluster up to Mpc scales. This is quite similar 
to what can be done with deep X-ray maps for interme- 
diate redshift clusters ( Boucail et al., 2000] ). We want to 
insist on the fact that the magnification bias effect is easy 
to detect from the observational point of view because it 
is less sensitive to seeing conditions or to geometrical dis- 
tortions of the instruments than are shear measurements 
( [Broadhurst et al., 1995| ; |Fort et al., 1997| ) . It only requires 
deeper observations, not necessarily in photometric condi- 
tions, provided one is able to reach the outer parts of the 
cluster to normalise the number counts. In order to im- 
prove the mass reconstruction, one may progress with the 
help of photometric redshifts, quite useful for faint objects, 



to constrain the redshift distribution of the sources ( Pcllc 
et al., 1999| ; Bolzonclla et al., 200C). This approach might 



uncertain, but a positive signal does not seem to extend 
much above ~ 1 ". More quantitatively, we can use the 
most recent measures of th e two point correlat ion func- 
tion from the HDF-South ( |Fynbo et al., 2000|) . At the 
magnitude limits used in our paper, this correlation func- 



also be useful in order to limit the influence of large scale 
over-density fluctuations when evaluating the magnifica- 
tion and the asymptotic limit of the depletion curves. This 
of course requires deep multi-color photometry, such as the 
one obtained on MS1008-1224. Extending this study up 
to large-scale structures (LSS) would probably be more 
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difficult to implement as compared to the search for cos- 



mic shear (Van Waerbeke et al., 200C), and the depletion 
effect should be restricted to cluster scales, at least with 
the simple method used in this paper. 

Finally, one of our initial prospects was to try to con- 
strain the background redshift distribution with multi- 
wavelength observations of the depletion effect. We have 
shown that this is quite a difficult task as the wavelength 
dependence of the depletion curves is a kind of second or- 
der effect. Nevertheless, this may be an interesting point 
to explore at other wavelengths, where the background 
redshift distribution is quite different than in the opti- 
cal. For example, dee p ISO observations i n the mid-IR of 
a few cluster lenses ( Altieri et al., 1999 ) may represent 
an extension of our analysis, as well as submm observa- 
tions with SCUBA, provided enough sources are detected 
behind the lenses for a statistical analysis. Another pos- 
sibility would be to address the question of the nature 
of the X-ray background sources and their redshift distri- 
bution through deep and high resolution observations of 
clusters with the new X-ray satellites Chandra and XMM 



(Rcfregier and Loeb, 1997) 
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